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Abstract 

We consider the measure problem in standard slow-roll inflationary models from the perspective 
of loop quantum cosmology (LQC). Following recent results by Ashtekar and Sloan, we study 
the probability of having enough e-foldings and focus on its dependence on the quantum gravity 
scale, including the transition of the theory to the limit where general relativity (GR) is recovered. 
Contrary to the standard expectation, the probability of having enough inflation, that is close to 
one in LQC, grows and tends to 1 as one approaches the GR limit. We study the origin of the 
tension between these results with those by Gibbons and Turok, and offer an explanation that 
brings these apparent contradictory results into a coherent picture. As we show, the conflicting 
results stem from different choices of initial conditions for the computation of probability. The 
singularity free scenario of loop quantum cosmology offers a natural choice of initial conditions, 
and suggests that enough inflation is generic. 
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I. INTRODUCTION 



The measure problem in cosmology has received some attention since it was suggested 
that one should weight, over the space of classical solutions to the equations of general 
relativity, those solutions that exhibit enough inflation to account for present observations 
PP. An early observation was that there exists a natural measure on the phase space of 
the theory with respect to which one should compute probabilities. Recently, Gibbons and 
Turok overcame some early difficulties in the total normalization and concluded that, for the 
simplest inflationary potentials in a FRW universe, the probability of inflation was greatly 
suppressed [2j. One potential difficulty with such calculations pertains to the choice of 
initial conditions. Since all solutions to the equations of motion are singular in the past 
(for expanding universes), one needs a prescription for selecting initial conditions for those 
solutions. In [2] such a prescription was put forward in terms of a 'constant density surface', 
roughly speaking, at the end of inflation. Another possibility is given by defining a 'cut-off', 
in the form of a constant density surface at, say, the Planck scale, as was early suggested in 

Yet another possibility is that a quantum theory of gravity might be able to provide 
such Planck surface in a natural way. Such is indeed the case of loop quantum cosmology 
[I] , a quantum framework closely related to loop quantum gravity [5] that has been able to 
achieve robust results regarding avoidance of big bang singularities P, [7J (See, for instance, 
[S] for a recent survey). In LQC, all trajectories undergo a bounce that replaces the initial 
singularity, attain a maximum critical density \7\. and preserve semiclassicality across the 
bounce pj, thanks to uniqueness results that warranty the consistency of the theory [1 Oj . 
Two key results in the measure problem have been obtained in LQC. First, it has been 
shown that one could account for the dynamics of the quantum universe by means of effective 
equations that capture the main quantum gravity effects and that reduce to the classical 
equations in the appropriate regime [HI [12]. This was used in [11] to show that, for several 
inflationary potentials, the characteristic 'attractor behavior' of inflationary dynamics [3|[T3l 
H3] is recovered in the low energy regime. Furthermore, Ashtekar and Sloan showed recently 
that the natural measure of [2] can be finitely implemented in LQC, and proposed a natural 
Planck scale surface on which to compute probabilities |15j . Surprisingly, the probability for 
having enough e-foldings was shown to be close to one, in contrast to the result of Gibbons 
and Turok that was done for classical GR 1 . 

In loop quantum cosmology, the underlying discreetness of the quantum geometry man- 
ifests itself via a dimension-full parameter A. In the LQC literature it is standard to choose 
the value of A such that the minimum quantum of area corresponds to that found in LQG 
[3 |6]. But, if one considers this as a free phenomenological parameter of the theory, it is 
natural to ask whether in the limit A — > 0, where the loop quantum geometric effects dis- 
appear, one can recover the standard Wheeler-DeWitt quantum cosmology. This has been 
answered with different levels of sophistication [HI CO HE]- The authors of [B] showed that the 
difference equation governing the LQC dynamics reduces to the differential WDW equation 
in the large volume limit. Later, in [7] and [18], the limit A — ?- was studied and it was 
shown that one does recover the standard WDW and the GR limit in some regime. In 



There have been several previous attempts to study the issue of inflation within LQC. In [IB] the natural 
measure of [2] was considered but the effects of the bounce and superinflation were ignored. In |17j the 
issue of the measure was not considered and only a small part of the parameter space was explored. 
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the case of effective classical equations, in this limit one recovers the equations of general 
relativity. 

The purpose of this article is to explore the relation between loop quantum cosmology 
and general relativity, as we take the limit A — > 0, regarding the measure problem in slow 
roll inflation. More precisely, we would like to understand the apparent tension between 
the results of Gibbons and Turok, with those of Ashtekar and Sloan. If one starts with the 
analysis of [15] , that was done for a fixed value of A (of the order of the Planck scale) , and 
one takes the limit A — > 0, one might expect to recover the results of Gibbons and Turok. As 
we shall show in detail this expectation is not realized. Indeed, quite the opposite occurs. 
As the value of the discreetness parameter is decreased, the probability of having enough 
inflation increases and approaches one in the limit. One would then be forced to conclude 
that in the general relativity limit of loop quantum cosmology, the probability of having 
enough inflation is (almost) one, in stark contrast with the analysis of Gibbons and Turok. 

What is then the source of this apparent tension? As we shall argue, the tension is 
resolved once one analyzes in detail the assumptions underlying both calculations. The 
difference turns out to be due to the initial conditions one imposes on the corresponding 
'constant density surface'. In the Gibbons and Turok analysis this is taken near the end 
of inflation, well below the Planck scale, whereas in the LQC calculation one is taking it 
at the scale set by the parameter A (which in the Ashtekar and Sloan analysis is close to 
the Planck scale). In the limit A — > the energy density at which the initial conditions are 
defined in LQC diverges, so one comes closer to the big bang singularity as one approaches 
the GR limit. It is this difference what accounts for the conflicting conclusions. 

The structure of the paper is as follows. In Sec. [II] we give a brief review of the effective 
description for loop quantum cosmology of a k=0 FRW cosmology with a scalar field. In 



Sec. Ill we present the calculation of the probability for having M e-foldings or more in 



LQC. We put special attention to the discreetness parameter of LQC and the limit when it 
vanishes. Next, we give an argument based on global properties of the dynamics and the 



Liouville measure to understand the results of both [15] and [2]. We end in Sec. IV with 
a discussion. Throughout the paper we use Planck units, where G=h=c=l, (rather than 
87rG =1, a convention sometimes used in cosmology). 



II. EFFECTIVE DYNAMICS IN LOOP QUANTUM COSMOLOGY 

Let us now give a brief review of the effective formalism in LQC. The effective Hamiltonian 
that one obtained from loop quantum cosmology for a k=0 FRW model is [T2] 

3 

W eS =-- — ^-v sin 2 X/3 + pv (1) 

where v is the volume and, on equations of motion, (3 = '-/H, where H is the Hubble 
parameter. From the previous Hamiltonian the effective Friedman equation becomes, 

sin 2 A/3 8n 

W = T" < 2 > 

or, equivalent ly 

ff 2 = (l - (3) 

3 V Pcrit / 
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FIG. 1: Three sets of trajectories are plotted for different values of the critical density. Note that 
near the origin, all trajectories approach the attractor. 

where the density is given by p = 4> 2 /2 + V(4>). Here p CTit = 3/(87T7 2 A 2 ), the critical density, 
is the density of the scalar field at the bounce. All trajectories undergo a bounce for which 
the density becomes exactly p C rit- I* 1 t ne low density regime, namely when A/3 < 1 or 
P "C Pcrit we approach classical general relativity. Note that the quantum geometry scale A 
sets the scale for the critical density. With the standard value taken in the LQC literature 
A = a/ 4v / 3tt7 ipi and the Barbero-Immirzi parameter 7 « 0.237 chosen to be compatible 
with the Hawking-Bekenstein entropy [19J, the critical density is p cr i t ~ 0.41pp/ [?]. (Recall 
that is the Planck units we are using £p t =Gh=l, and ppi = 1.) As we decrease the parameter 
A, the critical density increases, so the 'classical limit' is attained in the limit when the critical 
density diverges. 

The equation of motion for the scalar field yield the standard Klein Gordon equation 

is, 

4> + 3H<P + = (4) 

For the simplest potential, namely V = m 2 2 /2, we have solved the equations of motion 
for various values of the critical density and for convenience, plotted them in Fig. [TJ In the 
((p, <f>) plane, the surfaces of constant density are ellipsoids defined by p = <fi 2 /2 + m 2 (j) 2 /2. All 
trajectories approach the 'critical density surface', the ellipse bounding the phase diagram 
where the bounce occurs and touch it tangentially. Something that one might expect and 
that was checked in [11], is that near the origin of the plane, where the density is small 
compared to the critical one, the LQC trajectories and the classical one should coincide. 
This can be seen in Fig. [TJ As one decreases A the critical density increases and the maximum 
ellipse defined by <fi 2 B /2 + m 2 (f) 2 B /2 = p CTit becomes larger. The classical limit (GR) can be 
approached as A — » 0. One has to note however, that this limit is somewhat discontinuous 
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[7J [18], since all LQC trajectories bounce, for all values of A, while there is no bounce in 
GR. In this particular sense, the GR 'limit', and correspondingly the big bang, corresponds 
to an 'infinitely large ellipsoid', or the point at infinity in the (0,0) plane (See Fig. l2J). 



III. PROBABILITY FOR SLOW ROLL INFLATION IN LQC 

This section has two parts. In the first one, we calculate the probability for slow roll 
inflation in LQC and consider the limit when the discreetness parameter tends to zero. In 
the second part, we use qualitative aspects of the dynamics to gain a deeper understanding 
of the results. 



A. Probability 

Let us now evaluate the probability for inflation as done in |15j . keeping track of the 
dependence on A. Without loosing generality, for the remainder of the article we shall focus 
on the sector of the solution space for which is nonnegative. Then, the Liouville measure 
dp when pulled back to the surface with constant (3 or equivalently with constant p, has the 
from [T5] . 

dp = V8n~fv(p - 1/(0)) d0 dv (5) 
We we further choose, as in [15], the surface of constant (3 (and p) at the bounce, we get 

dp = — — vbV 1 ~ f b d(t) B dv B (6) 

where 0# is the value of scalar field at the bounce, v B is the volume of the universe at the 
bounce and F B = V \4> B ) / Pent- This is the measure that will be used for computing the 
probability of having M or more e-foldings. 

The number of e-foldings during inflation, Af, can be written as 

M= Hdt= / -^d0 (7) 

Jt J<t>o 

where t , O , t en d, and en d are the time and value of the scalar field at the onset and at the 
end of inflation, respectively. We can use the slow roll conditions, V(4>) ^> 0/2 and V$ ^> 0, 
together with Eq. Q to approximate J\f, 

•Af « - ^ # = 2n ( 1 - ^fiA (<Pl ~ 0L) (8) 



V, 



max 



where max is the maximum value the scalar field can attain and is given by max — 
^2p crit /m. For large values of the value of the scalar field at the end of inflation is 
much smaller than its value at the onset of inflation. Thus, for large (but finite) M we can 
neglect some terms and get, 

■A^ 2tt (l - -JU 0* (9) 
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It should be noted that this is a slight overestimation of the value of M but this does not 
constitute a problem for our analysis. From this last equation, we can find the value 0^ of 
the scalar field at the onset of inflation, for a given value of M as 



^3(1 - 8A/yA 2 m 2 /3) 



¥L = ±^ — y - 7 ^~ — (io) 

Vo7T7Am 

In the GR limit, that is, in the A t— > limit, we expect that 0^j_ be equal to ±a/A/"/27t. 

Let us now see how we can find 4> B which is the value of the scalar field at the bounce 
that evolves under the dynamics to <p ± as the starting point of inflation. According to 
Eq.(4|, if at the bounce 4>° B > then <f>B < (and > 0). Similarly, if <p° B < then 0# > 
(and <pB > 0). In the second case, after some time, becomes zero and after that it will be 
negative, but near the onset of inflation it becomes zero again. Near the start of inflation 
at the time for which = 0, the value of the scalar field is larger than <ft° B . After that, 
becomes negative and the value of the scalar field starts to decrease but very soon after 
= 0, the inflationary era starts and the scalar field at the onset of inflation remains larger 
than the value of the scalar field at the bounce (0^ < O ±)- 

Furthermore, because of the uniqueness of the solutions, O is a monotonic function of 
4>b and since O is always greater than 0#, then it is an increasing function of 0s . 

Given this, we can write the probability of having inflation with M e-folding or more as 
the quotient of the volume on the space of solutions occupied by solutions with Af or more 
e-foldings divided by the total volume. Since the measure does not depend on volume v 
and the range of this coordinate is infinite, both terms are unbounded. However, we can 
very easily get rid of these spurious infinities by an appropriate renormalization (or gauge 
fixing [20]). One possibility is to restrict the domain of the volume integral to the interval 
v G (1,2) (in Planck units). With this choice, the volume integrals in the quotient cancel 
each other and we get, 



P 



/-1L V 1 - g + Jgr VT^F-B d0B 



fr v / r^d0 B 

where 0^ and 0^ are the minimum and maximum value of at the bounce that cause 
inflation with M e-folding respectively and max is the maximum value of 0s and is equal 
to 3/27Ami/7r. Then 



arcsin(2 7 Am0^ v ^F/3) - arcsin(2 7 Am0f 0r/3) - 2 7 Am v ^F/3(0^ - 0^) 
N ~ V 2(tt/2 - 1) {U) 



We have plotted in Fig. [T]the dynamical trajectories for three values of A. As one can 
see, when A becomes small, the trajectories (for finite values) are almost parallel. Then, in 
the limit A — > we can approximate — by 0^j_ — <p^_ and since 0^_ are finite, then the 



difference between 0^ and 0^" is finite. From the above discussion and Eq. ( 12 ) , for a finite M 
we see that the probability is a decreasing function of A (27Am- v / / 7r0f,/3 and 27Am v /7r0 a /3 < 
1) and when A goes to zero we have that arcsin(27Am0^ / ' \A"/3) ~ 2^j\m(j) J ^ (and 
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equivalently for (p^) and therefore the probability in Eq.(12) goes to 1. This is the first 
result of this paper. 

Let us now understand qualitatively why the probability increases as the LQC parameter 
decreases. As the analysis here presented and that of [TH] shows, for a given value of p crit , 
there is an interval (0^,0^0, in the 'kinematically dominated regime' (where the energy 
density at the bounce is mainly due to the kinetic energy), where there are not enough 
e-foldings. This interval, as we have estimated before, depends on A. In Fig. [2] we have 
plotted, for three values of p cr it, the 'critical trajectories' for which the transition occurs. 
That is, these trajectories have an almost identical behavior at small densities, so they 
inflate in the same fashion, and touch the bounce surface at the points <ft^ and ' . If we 
now follow them to higher densities 'back in time', what one sees from the graph is that 
as A decreases, and p crit increases, the intersection points tend to the axis. The relative 
size of the interval (0^, 0^) in the total allowed interval (— max , 0max) also goes to zero as 
p cr it — > oo 2 . Since the integrand does not diverge, this already implies that the quotient 
vanishes and the probability goes to 1. 

Note also that this result is independent of the precise value of M (as long as it is large 
enough for our approximation to be valid). Does this mean that we can take the limit 
M — > oo and also have probability one? In order to answer this one should exercise some 
care. For any finite value of A/", the probability will tend to one as we make / smaller for 
two reasons. The first one is that the dynamics of the effective equations is such that those 
trajectories that do not have enough inflation get 'funneled', for large enough values of the 
critical density, into the interval (0^", ) that remains bounded, while the total interval 
for grows with p cr it- This only happens because we are taking the bounce surface as the 
reference surface where the probability is computed. Furthermore, the measure is such that 
relative volume we associate to those trajectories is very small and becomes zero in the / — > 
limit. This does not mean that we can fix I and, say, take the limit N — > oo. 

A final remark is in order. In our analysis, as plotted in Fig. |2j the criteria for how much 
inflation there is coincides with that of [15] . That is, we start with the critical density of 
LQC (of the order of the Planck density), which gives an initial condition from which to 
measure e-foldings, and find those trajectories -in theories with a different I- for which the 
dynamics at low densities coincide, where inflation actually occurs. This is also in the spirit 
of [3], which suggested to take initial conditions at the Planck scale. Our strategy has to 
be contrasted with a possible alternative that involves going closer to the big bang, as we 
decrease /, and use that as initial condition in the e-folding counting. The problem with this 
choice is that, as one approaches the big bang that has zero volume, the number of e-foldings 
diverges for all trajectories, so even the question of which trajectories have enough inflation 
becomes meaningless, since every trajectory would have an infinite number of e-foldings. 



B. Comparison 

Let us now come to the question of how we can reconcile the results of Gibbons and 
Turok [2] on the one side and those of Ashtekar and Sloan [15] on the other. The first 
possible objection is: How can we compare two results that are taken on two different 



2 Recall that </> max is obtained from the value of the potential at the bounce p cr it = V(<fr max ). In our case 
<^max = \/2 Pcrit /m, so the interval in which the scalar field can take values also diverges as p cr ;t — > oo. 
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FIG. 2: Here we are plotting trajectories for three different values of the critical density p CT it- In 
each case, we have a boundary of the trajectories in the (</>, cj>) plane, corresponding to the bounce, 
that are depicted as ellipsoids. The smallest ellipsoid can be taken as the LQC one, and the larger 
ones are closer to the GR limit. The 'critical' curves that separate the region of enough e-foldings, 
as determined by the LQC scale, are then plotted for the three different values of the critical 
density p cr it- One can see that, as the critical density increases, the intersection with the ellipsoid 
of critical density comes closer to the (ft axis. 

theories, GR on one side and LQC on the other? As we have seen before, one can in 
fact approximate very well the low density GR trajectories by (low density) LQC effective 
trajectories. Thus, the region of interest in the Gibbons and Turok analysis, per/Pent "C 1, 
which is for trajectories near the end of inflation (and therefore, around the constant density 
surface in our Figure [2j, one can take the LQC effective trajectories without any problems as 
a very good approximation to the GR dynamics. This allows us to 'embed' the low density 
GR dynamics in the effective LQC description with very good accuracy. 

With this assumption, we can now compare the two result within the effective LQC 
description. We have two constant density surfaces, as depicted in Fig. [3] The external 
ellipsoid corresponds of course to the critical density p crit at scale A, while the small one 
corresponds to the density pgt as chosen by Gibbons and Turok 3 . One puzzling fact about 
the huge discrepancy in results is that both analysis use the natural Liouville measure (prop- 
erly normalized) to compute the probability on constant density surfaces. One important 
property of the Liouville measure is that it is invariant under the dynamical evolution. So, 
how come we arrive to two very different conclusions? 



3 The figure is not to scale, since we are asking that pgt "C Pcrit- The relative densities in the figure were 
taken to illustrate our point. 
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There are two key observations to understand this apparent tension. The first one pertains 
to the question of whether the time evolution invariance of the Liouville measure implies that 
the probability is also invariant. On a first view, one might imagine that the probability 
has to be invariant since one is just measuring the relative phase space volume of those 
trajectories with M e-foldings or more, relative to the total volume in phase space. Now, 
the technical step that allowed to normalize the phase space volume (the total phase space 
volume is infinite) in [2] and [15] was to realize that there is an invariance in the space of 
classical solutions by rescaling the physical volume. This invariance has its origin in the 
fact that, instead of describing the whole universe, one has to restrict attention to a fiducial 
region R in space (the spatial volume of the whole universe in k=0 FRW is infinite, so one 
needs to consider a region with a finite volume). Since this choice is arbitrary, one can in 
principle chose a smaller /larger region for which we assign a smaller/larger volume, but the 
physics should be unchanged. When one takes care of this ambiguity, either by taking an 
appropriately chosen 'interval in u' as we done in the previous part, or by an appropriate 
gauge fixing [151 ED] , one still has to be careful about the possible change in physical volume 
during the dynamical evolution that would also induce a change in relative volume in phase 
space. 

Let us see how this comes about. Invariance of the Liouville measure means that the 
volume in phase space is preserved. Let focus our attention in the quadrant in the space 
of solutions, with coordinates (v B ,<p B ), defined by 1 < v B < 2 and — max < <Pb < 0max ? 
and follow it through its dynamical evolution. If we now take another 'gauge fixing' at a 
lower energy density, say p 1 (See Fig. [3]), we immediately notice that the range in <p is much 
smaller. Since the total volume of the quadrant we are following has to be the same (due 
to the dynamical invariance of the Liouville measure), the range in v has to increase, as it 
indeed does, since most solutions inflate. The crucial point here is to realize that the change 
in volume Av = v i — v b from the bounce to the pi surface depends on the value of (p. Thus, 
the lines, say, vb(4>b) = 1 at the bounce gets mapped, in general to a curve t>i(0i) that is no 
longer constant as a function of <pi. That is, each solution has a different change in physical 
volume depending on the value of the scalar field at the bounce. But, if one is only keeping 
track of the change in the 'phase space coordinate' <fi when computing the probability, then 
the relative volume in phase space, as measured by only 0, can indeed change. Since this 
is precisely what one means by probability in the analysis of [2j [15] and here, we are led to 
conclude that the probability indeed depends on the surface on which it is computed. Since 
this argument did not use any particular detail of the LQC dynamics, this ambiguity in the 
probability depending on the choice of constant density surface is also present in general 
relativity. Let us now see what further assumption are made in both calculations. 

The second observation is the following. When computing the probability of having M 
e-foldings, one has to assume an a-priory probability distribution P(0, v) of the classical tra- 
jectories, and then integrate this probability distribution with respect to the corresponding 
measure. In [15], the authors consider the most natural choices, namely, the probability is 
computed on the critical density surface (i.e., the bounce) using, as the integration measure, 
the Liouville measure. By invoking Laplace's 'principle of indifference' as in [2], they con- 
sider a uniform distribution on the space of trajectories (labeled by (<fi,v)) and performed 
an appropriate gauge fixing with respect to the volume rescaling freedom available, in the 
same spirit we have done here. We have illustrated this scenario in Fig. [3j where we plotted 
trajectories uniformly distributed in <ft along the critical density surface. If we now follow 
these trajectories along the dynamical evolution we notice that, when they intersect the 
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FIG. 3: For a fixed value of p C rit, we plot the exterior, critical density surface and a surface of 
constant density pgt *C p C rit (not drawn to scale, of course) on the (0, 0) plane. Trajectories with 
a uniform distribution at the LQC bounce ellipsoid are plotted. Note that trajectories for which 
there is enough inflation get funneled into a small region in the smaller /9qt ellipse. Near this 
surface, the GR and LQC dynamics almost coincide 

P = Pgt surface, they are no longer uniformly distributed. Quite the opposite. Due to the 
global properties of the dynamics, the trajectories are funneled into the 'attractor' on the 
plane (0,0) and, therefore, effectively acquire a new 'probability distribution' P(0) on the 
P = Pgt surface 4 . 

If we were to compute the probability of inflation on the Gibbons and Turok surface, 
but weighted with the induced distribution P(0), we would of course get the same result of 
Ashtekar and Sloan, given our previous discussion. What Gibbons and Turok did instead 
was to assume a uniform distribution P(0) on the p = pgt surface. With respect to the 
uniform distribution, the phase space volume of inflating solutions is very small and the 
probability is therefore, very close to zero. Had we chosen to compute the probability on a 
surface with even lower density, the result would even be smaller. This constitutes the main 
difference between both calculations 5 . 



4 We can view this induced probability distribution as a way of keeping track of the relative change in 
phase space volume do to the dynamics that induces a differentiated change in physical volume v for 
different trajectories. Namely, the distribution is given by P(0) = vgt(4>), in terms of the volume i>gt, 
as a function of <j> on the Gibbons- Turok constant density surface. 

5 Note that Linde had pointed out that the assumption of Gibbons and Turok to take uniform initial 
conditions at the end of inflation might be the source of their negative result P3] , but the mechanism he 
outlined is different from ours, since he was not taking the Liouville measure into account. 
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Furthermore, we can now understand why the probability found by Gibbons and Turok 
is so small. If we look at the region of Fig. [2] for which there is not enough inflation (in 
between the critical curves) on the LQC critical density surface, and follow those trajectories 
as in Fig. [3j we see that those trajectories occupy now a much larger region on the p = pgt 
ellipsoid. In other words, the trajectories for which there is enough inflation get funneled 
into a small region in the p = pgt surface that, when integrated with respect to the uniform 
distribution P(0) of [2], give a very small contribution to the probability. One could also 
consider the opposite situation in which one starts with an uniform distribution on the GT 
surface and 'evolve back' in time to the bounce surface. In that case, the dynamics will 
'expel' the trajectories in such a way that the probability distribution P'(0) induced on the 
bounce surface is concentrated on the region where there is not enough inflation. If one 
integrates that probability distribution with respect to the Liouville measure the resulting 
probability is very close to zero, as found by Gibbons and Turok in [2|. 6 

IV. DISCUSSION 

Let us summarize our results. We have reanalyzed the treatment of the simplest infla- 
tionary model from the perspective of loop quantum cosmology. By using effective equation 
we studied the structure of the space of classical solutions with the aim of answering the 
question: How probable is it to achieve enough e-foldings? In particular we have consid- 
ered this question keeping the discreetness parameter of loop quantum cosmology as a free 
parameter. When the parameter vanishes, one expects the dynamics to reduce to the stan- 
dard, general relativity behavior. The first result is that, as previously shown in [15], the 
probability for enough inflation is very close to one when the discreetness parameter I is of 
the order of the Planck scale. We then considered the dependence of the probability as one 
decreases the parameter and it approaches the general relativity limit. As we have shown, 
the probability increases and approaches one as one reaches the limit. Next, we studied the 
global properties of the system to understand the underlying reason for the discrepancy of 
these results and those of [2] in which the probability of enough inflation was computed to 
be close to zero, within general relativity. What we found is that this discrepancy is due 
to the differences in the underlying assumptions in both calculations. As it turns out the 
probability as computed in both [2] and [15] depends very strongly on the constant density 
surface where it is calculated. While Ashtekar and Sloan assume a uniform distribution of 
classical trajectories at the naturally defined surface available due to the universal existence 
of the bounce, Gibbons and Turok take it at an arbitrarily defined surface at the end of 
inflation. Given the large difference in scales involved and due to the global properties of 
the dynamics and the probability measure, these two assumptions have strikingly different 
consequences. During the evolution from the bounce to the Gibbons- Turok scale, most of the 
trajectories that undergo enough inflation -contributing significantly to the probability- get 
funneled into a small region at the later scale, that has a correspondingly small contribution 
to the probability. This is the origin of the apparent tension. 

We have thus found two very different results even for GR. On the one hand the Gibbons- 



6 One should note that, as previously discussed, in order to make the distinction of which trajectories have 
enough e-foldings, one has to introduce a cut-off for the initial condition. Here we have adopted the LQC 
scale as a natural unambiguous choice. 
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Turok result involves several, somewhat ad-hoc, assumptions given that there is no preferred 
choice of scale on GR. On the other hand, there is the limit of LQC when the discreetness 
parameter vanishes. In this later case we have, for each scale, calculations based on unam- 
biguous and natural choices that provide a well defined result, even when the GR limit of 
LQC is non-smooth. Thus even when in the / — > limit one is approaching arbitrarily close 
to the singularity, the probability of having M e- foldings -as measured from the Planck scale 
down- can be given some meaning. As we have seen, the result that in the GR, I — > limit, 
the probability goes to one for any finite value of J\f, seems to be generic. Whether this 
result is physically meaningful is, however, a completely different issue. In the situation in 
which A is taken well below the Planck scale, we are implicitly assuming that the classical 
equations are still valid. This is perhaps, too strong an assumption. One generically expects 
quantum effects to dominate near and below the Planck scale. This is precisely what LQC 
provides for us via its effective equations. 
Let us end with a series of remarks. 

1. Given that these two results are based on assumptions that yield completely oppo- 
site predictions, one might then ask what is the physically reasonable assumption to 
make? How can we justify one choice over the other? Is there a 'canonical' choice of 
initial condition? In loop quantum cosmology we know that the bounce is generic for 
inflationary potentials |21| . and the effective equations are a very good approximation 
to the dynamics of semiclassical states [12J. Since every such effective trajectory goes 
through a bounce, selecting the surface of constant density, at the bounce, seems a 
rather natural choice. One should emphasize then that there does not exist a similar 
surface that is preferred in the classical GR case. As we have seen, the probability 
does depend in a rather dramatic way on the choice of such surface. Without any 
extra input, the LQC choice seems to be the most natural. 

2. Even if one does not regard loop quantum cosmology as a fundamental theory, one 
can still view its effective dynamics and choice of surface as in [TBj as a procedure 
to regulate the classical calculation. The bounce provides then the preferred 'cut-off' 
surface envisioned by the authors of [3], but in an unambiguous fashion. From this 
perspective, what is amazing is that one can remove the regulator and obtain a finite 
answer. Furthermore, this 'canonical answer' indicates that inflation with enough 
e-foldings is generic even in this particular way of approaching the GR limit. 

3. As we have seen from our analysis here, the reason for the LQC result stems from the 
choice of surface where to compute the probability and not directly from any effect 
from the quantum geometry underlying LQC. In fact, as we remove the parameter 
encoding the quantum geometric effects, the probability of inflation increases 7 . 

4. In our analysis we have used qualitative aspects of the global dynamics of the sys- 
tem. Therefore, our arguments and conclusions are insensitive to changes in the free 
parameters of the model. That is, our result are rather robust. 

5. One should keep in mind that these results are purely classical. It is to be expected that 
quantum effects might provide a more realistic distribution on the space of classical 



7 This is consistent with early calculations on LQC where the integration was performed at low densities 

EH- 
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trajectories. Some proposals have been put forward, but in the context of particular 
states and only for those trajectories satisfying WKB conditions [22]. One could 
imagine that semiclassical states in LQC might provide an improved distribution from 
which one might get a 'quantum corrected' estimation of the probability for enough 
inflation. This matter should certainly be studied in detail. 
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